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During the core bounce of a supernova collapse resonant active-to-active (y a — » h> a ), 
as well as active-to-sterile (u a — * u s ) neutrino (y) oscillations can take place. Besides, 
over this phase weak magnetism increases antincutrino (i>) mean free paths, and thus 
its luminosity. Because the oscillation feeds mass-energy into the target v species, the 
large mass-squared difference between species (u a — » v s ) implies a huge amount of power 
to be given off as gravitational waves (£q"Ws ~ 10 49 erg s — 1 ), due to anisotropic but 
coherent v flow over the oscillation length. This anisotropy in the i^-flux is driven by both 
1 the universal spin-rotation and the spin-magnetic coupling. The new spacetime strain 

$_H ' estimated this way is still several orders of magnitude larger than those from v diffusion 

(convection and cooling) or quadrupole moments of the neutron star matter. This new 
feature turns these bursts the more promising supernova gravitational-wave signal that 
may be detected by observatories as LIGO, VIRGO, etc., for distances far out to the 
VIRGO cluster of galaxies. 
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1. Generating GWs During v Oscillations 

To start with, let us recall how the production of GWs during v oscillations proceeds 
by considering the case of oscillations between active and sterile neutrinos in the 
supernova core. The essential point here is that oscillations into sterile neutrinos 
change dramatically the energy and momentum (linear and angular) configuration of 
the proto-neutron star (PNS). In particular, flavor conversions into sterile neutrinos 
drive a large mass and energy loss from the PNS because once they are produced 
they freely escape from the star. This means that oscillations into sterile, in dense 
matter, take place over longer oscillation lengths, compared to v a — * v a , and the 
sterile encounter infinite mean free paths thereafter. Physically, the potential, V s (x), 
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for sterile neutrinos in dense matter is zero, and their probability of reconversion, 
still inside the star, into active species is quite small. This outflow translates into a 
noticeable modification of the PNS mass and energy quadrupole 

Since most sterile neutrinos escape along the directions defined by the dipolc 
field and angular momentum vectors (see Fig. 2), then the v outflow is at least 
quadrupolar in nature. This produces a super strong gravitational-wave burst once 
the flavor conversions take place, the energy of which stems from the energy and 
momentum of the total number of neutrinos participating in the oscillation process a . 

The remaining configuration of the star must also reflect this loss, i.e., the active 
vs were coupled to the neutron matter. Hence its own matter and energy distribution 
becomes also quadrupolar. Because this quadrupole configuration (the matter and 
energy still trapped inside the just-born neutron star) keeps changing over the 
timescale for which most of the oscillations take place, then GWs must also be 
emitted from the star over that transient. At the end, the probability of conversion 
and the z^-flux anisotropy parameter (a, see below) determine both how much energy 
partakes in the process and the degree of asymmetry during the emission. 

The case for oscillations among active species is a bit different, the key feature 
being that mass and energy is coherently relocated (see Ref. 1 ) from one region to 
another inside the PNS, especially because of the weak magnetism of antineutrinos 
that allows them to have larger mean free paths 2 . In addition, oscillations of electron 
neutrinos into muon or tauon neutrinos leave these last species outside their own 
neutrino spheres, and hence they are in principle free to stream away. This must 
also generate GWs during that sort of flavor conversions, although their specific 
strength (strain) must be a bit lower compared to conversions into sterile neutrinos 
wherein almost all the v species may participate, and the large Am 2 in the process. 
These neutral-current interacting v species must be the very first constituents of 
the v burst from any supernova since most ^ e s are essentially trapped. 

Since the sterile neutrinos escape the core over a timescale of a few milliseconds 
(ms), the number of neutrinos escaping and their angular distribution is sensitive 
to the instantaneous distribution of neutrino production sites. Since thermalization 
cannot occur over such short times (AT* hermal ~ 0.5 s), and since the neutrino 
production rate is sensitive to the local temperature at the production site, the 
inhomogeneities during the collapse phase get reflected in the inhomogeneities in the 
escaping neutrino fluxes and their distributions. Because of both the v spin-i? and 
v spin- J coupling the asymmetries in these distributions can give rise to quadrupole 
moments, which must generate gravitational waves 3,4 , and dipole moments which 
can explain the origin of pulsar kicks 5 ' 6 . 

Fixed by the probability of oscillation, P Va ^ Va , the fraction of neutrinos that can 
escape in the first few milliseconds is, however, small. Firstly, the neutrinos have 
to be produced roughly within one mean free path from their resonance surface. 

a The attentive reader must regard that neutrinos carry away almost all of the binding energy of 
the just-born neutron star, i.e., AE„ ~ 3 X 10 53 erg. 
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Secondly, since in the case of v a — > v s oscillations the v s is the heaviest neutrino 
species, the sign of the effective potential A(x) and the resonance condition indicates 
that only v e s and the antineutrinos D u and v r can undergo resonant conversions. 

2. Anisotropic Neutrino Outflow: origin and computation 

To provide a physical foundation for the procedure introduced here to determine 
the neutrino asymmetry parameter, a, which measures how large is the deviation of 
the v-fhix from a spherical one, we recall next two fundamental effects that run into 
action once a PNS is forming after the SN collapse. We stress in passing that other 
physical process such as convection, thermalization, etc., are not relevant over the 
timescale under consideration: AT osc ~ 3 — 10 ms after the SN core bounce, and 
therefore do not modify in a sensitive manner the picture described below. Indeed, if 
v thermalization, for instance, already took place, then the oscillations are severely 
precluded since oscillations benefit of the existence of energy, matter density and 
entropy gradients inside the PNS 7 ' 8 , which are "erased" once thermalization onsets. 

As discussed by Mosquera Cuesta and Fiuza 1 (see also the detailed discussion 
by Loveridge 6 about the v distribution function inside the PNS, an asymmetric, 
off-centered rotating v beam) both the i/-spin coupling to rotation and the f-spin 
coupling to the magnetic (B) field drives an effective momentum, and thus energy 
flux, asymmetry along the angular momentum (J) and -B-field directions. The com- 
bined action on the escaping us of a rotating background spacetime plus magnetic 
field makes their ^-sphere a decidedly distorted surface, with at least a quadrupolar 
distribution. This is the source of the strong GWs burst in this mechanism when 
the oscillations ensue. 

Based on the concomitant action of both the effects: j/-spin coupling to both 
the magnetic field and rotation, described previously, one can determine the v flow 
anisotropy in a novel, self-consistent fashion by defining a as the ratio of the total 
volume filled by the distorted ^-spheres to that of the proto- neutron star (PNS). The 
z^ e -sphere radius of a non-magnetic non-rotating star is obtained from the condition: 



where T„ e is the optical depth, and K, Vc the scattering opacity for electron neutrinos, 
and p the matter density. 

One can take hereafter 9 i?p c NS = R v ~ 35 km, which is of the order of mag- 
nitude of the oscillation length A„ of a typical supernova v\ which is constituent 
as well of the atmospheric z/s, for which Am 2 <~ 10~ 3 eV 2 has been estimated by 
Superkamiokande v detector 10 



Therefore, resonant conversions between active species may take place at the posi- 
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tion r from the center defined by 

r = R Ve + 5o cos <f> , 
with (j) the angle between the ^-spin and B, i.e., 

COS = ; s = -jh 73 x ~ 1 - 10 km , 

fc 2ir I (dN e /dr) 

for <~ 10 14 ~ 15 G, respectively. Here e, iV e = Y e N n and represent the electron 
charge, density and chemical potential, respectively. 

This defines in Fig. 2 an ellipsoidal figure of equilibrium with semi axes a = 
R v + 5q\ and b = R v , and volume (after rotating around B): 

Veiups. = ^nRl(Ru + <5 ) . 

Meanwhile, the ^-spin coupling to rotation generates an asymmetric lcmniscate- 
like plane curve (see Fig. 2) r = R v (1 ± L a cos 8) 1 , which upon a 2ir rotation 
around the star angular momentum axis generates a volume: 

Viemnisc. = \{1R v f x 2ny + ^{R v f x 2iry , 

where the quantity y (and x) is defined as the location of the centroid of one of the 
lobes of that plane figure with respect to its coordinate center (x,y). After a long, 
but straightforward, calculation one obtains y = r2 "^ 4 ^ R v , x — r 2^y 4 ) Rv ■ 
Thence, for a PNS 9 with parameters cited in this article one obtains 11 

Van,,. + Vi 

( in a /sc.. n 1 1 n m /-. \ 

a min = 1 -— 0.11-0.01. (1) 

VPNS 

As such, this is essentially a new result of this paper. The attentive reader 
must notice in passing that the definition in Eq.(l) does take into account all of 
the physics of the neutrino oscillations: luminosity, density gradients and angular 
propagation, since Fig. 2 does gather the relevant information regarding the spatial 
configuration of the v luminosity in as much as is done in the standard definition 
of a 9 ' 12 

Indeed, one can get the "feel" of the relationship between these two definitions 
by noticing that the quantity L v (Q',t) in the integrand of Eq.(2) can be expressed as 
L v (Q',t) = Lu(t)F(n'), where the function F(Sl') contains now all the information 
regarding the angular distribution of the neutrino emission. Hence L v (t) can be 
factorized out of the integral and dropped from Eq.(2). This converts Eq.(2) in 
a relationship among (solid) angular quantities, which clearly can be reduced to 
a volumetric one, similar to the one introduced in Eq.(l), upon a transformation 
using the definition of solid angle in the form of Lambert's law: dfl = dA ^° s 8 , and 
applied to the sphere representing the PNS. Here 9, measured from a coordinate 
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system centered at the PNS 12 , plays the role of the angle between the direction 
towards the observer and the direction ft' of the radiation emission in Eq.(24) of 
Ref. 12 . Therefore, the novel result here presented is physically consistent with the 
standard one for the asymmetry parameter a. 



Pulsar 




Fig. 1. A schematic representation of the !/-fiux distribution (hatched regions) in a nascent pulsar. 
The at least quadrupolar distribution is evident and stems from the neutrino spin-magnetic and 
spin-rotation couplings. 



3. Enlarged u and GWs Luminosity From Oscillations in Dense 
Matter 

Numerical simulations 12 showed that, in general, the fraction of the total binding en- 
ergy emitted as GWs by pure v convection 3,11,12 is: E^ w ~ [10~ 10 -10~ 13 ] M Q c 2 , for 
a v total luminosity: L v <~ 10 53 erg s _1 . In the production of GWs via v oscillations 3 
[y a <-» v a or v a <-» v s ) there exists two main reasons for expecting a major enhance- 
ment in the GWs luminosity during the transition: a) the conversion itself, which 
makes the overall luminosity (L Vx ) of a given v x species to grow by a large factor: 
AL;/ X < 10%L^°^ a ^. The enhancement stems from the mass-energy being given to, 
or drained from, the new v species into which oscillations take place. This augment 
gets reflected in the species mass-squared difference Am 2 , and their relative abun- 
dances: one species is number-depleted while the other gets its number enhanced. 
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But, even if the energy increase, or give up, is small, b) the abrupt (resonant) con- 
version over the transition time ATbsc = T -- A ° sc ~ [10~ 4 — 10 _3 1 s, also magnifies 

Vv-Diff. _ " 

transiently L Vx . Here Aosc defines the oscillation length, and V v -Diff. ~ 10 9 cms -1 
the convective v diffusion velocity 12 . 

If flavor transitions can indeed take place during supernovae (SNe) core collapse 
and bounce, then they must leave some imprints in the SNe neutrino spectrum. Main 
observational consequences of neutrino conversions inside SNe include the partial or 
total disappearance of the neutronization peak; the moment at which most u e s are 
produced, the interchange of the original spectrum and the appearance of a hard 
p e spectrum, together with distortions of the v e energy spectrum and alterations of 
the v e spectrum 13 . 

3.1. Resonance, Adiabaticity and the Role of Weak Magnetism 

v oscillations in vacuum produce no GWs 3,4 . In the case of active-to-active v os- 
cillations (essentially the same physical argument holds also for active-to-sterile v 
oscillations), the main reason for this negative result is that this class of conversions 
do not increase in a significant figure the total number of particles escaping from 
the proto-ncutron star. In dense matter, the process generates no GWs since the 
oscillations develop with the neutrinos having very short mean free paths, so that 
they motion can be envisioned as a standard diffusion process. 

However, because of the active antineutrino species weak magnetism their effec- 
tive luminosity can be enlarged as much as 15% compared to the typical one they 
achieve when this effect is not taken into consideration during their propagation in 
dense matter 2 . This result can be interpreted by stating that the number of oscil- 
lating (and potentially escaping) antineutrinos may be augmented by a large factor 
because now they do encounter longer mean free paths. Below we take advantage 
of this peculiar behavior of the v outflow in supernovae in computing the overall 
probability of transition between active species. 

3.2. Active- To- Active v Oscillations 

The neutrino oscillation pattern in vacuum can get noticeably modified by the 
passage of neutrinos through matter because of the effect of coherent forward 
scattering 14 ' 15 . Therefore, interaction with matter may help in allowing more ^s 
to escape if resonant conversions into active 3 and/or sterile i^s 3 occur inside the 
^-sphere of the active vs. The transition probability 7 for conversions among active 
species thus reads 

p »a^va (l^i - £o|) = 1 - ^ [1 + cos26> m (x ) cos26» m (xi)] , (3) 

where xq and x\ correspond to the production and detection sites, respectively, and 
9 m the mixing angle in matter. 

In order to produce an effect, neutrinos must be able to escape the core without 
thermalizing with the stellar material. For active neutrino species of energies w 10 
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MeV, this is not possible as long as the matter density is > 10 10 g cm' 3 . Since 
the production rate of neutrinos is a steeply increasing function of matter density 
(production rate oc p n , where p is the matter density and n > 1), the overwhelming 
majority of the neutrinos of all species produced are trapped. This way, there seems 
to be no contribution to the GWs amplitude for neutrino conversions taking place 
within the active neutrino flavors. Nonetheless, if there were indeed weak magnetism 
effects 2 , one can rethink of conversions within active species. In the new result on 
weak magnetism for antineutrinos in core collapse supernovae, the antineutrinos 
(9 x s) luminosity could be noticeably increased because of their longer mean free 
paths 2 , and this means that the total energy flux can be augmented in ~ 15% for 
^s of temperature <~ 10 MeV. One can verify that longer mean free paths allows for 
a larger oscillation probability, and hence the contribution to the generation of GWs 
during flavor conversions within active species becomes non-negligible compared to 
the earlier case 3 where not weak magnetism effects were taken into account. 

For active-to-active v flavor conversions, for instance: v e — ► fcy,, ^v; as implied by 
SNO results, the resonance must take place at a distance ir ros from the PNS center 
and amid the active ^-spheres, whenever the resonance condition is satisfied 

4At7j cos 29 m = 2V2 G F N e (x rcs ) k Ve ; 

notice that we neglected the magnetic field contribution. Here k Ve — E Ve /c is the 
v e momentum, and N e (x Tes ) — N e - — N e + <~ io 39 ~ 40 cm~ 3 is the electron number 
density. Thus the right-hand part of this equation reduces to 

E u . 



2V2G F E Ve N e (x ies ) = 6.9 x 10 4 eV 2 



10 n g cm - 



10 MeV 



(4) 



For densities of order p x=Xrcs ~ 10 10 g cm -3 , i.e., recalling that v e s are 
produced at the PNS outermost regions 3 where the electron to baryon ratio is 
Y e ~ 0.1, the resonance condition if satisfied for Am 2 ~ 10 4 cV 2 , which im- 
plies a neutrino mass of about m„ ~ 10 2 eV. Neutrino flavor conversions in 
the resonance region can be strong if the adiabaticity condition is fulfilled 3 , i.e., 

whenever 7 ^ ^26™ ("if ) » 1 , where x ICS is the position of the reso- 

nance layer. Recalling that the typical scale of density variations in the PNS core 
is /ijv e ~ (dNe/dr) -1 <~ 6 km, this adiabatic behavior could be achieved as far as 
the density and magnetic field remain constant over the oscillation length 

f I dpV 1 / 1 Am 2 , „ A" 1 1 cm 

A osc =(-/■) ~ — ^r- sin 2^ ro ~ . . , 5 

of order /ijv e , which can be satisfied for Am 2 ~ 10 4 eV 2 as long as 19 sin 2 29 m > 10 -8 . 
Although these v oscillations could be adiabatic for a wide range of mixing angles, 
and thus a large number of ^s could actually oscillate, a v x mass such as this is 
incompatible with both viable solutions to the Solar Neutrino Problem (SNP) and 
the most recent cosmological constraints on the total mass of all stable neutrino 
species that could have left their imprint in the Cosmic Microwave Background 
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Radiation (CMBR), as inferred from the observations performed by the satellite 
WMAP: m v ~ 1 eV. Therefore, we dismiss this possibility since there appears to 
be no evidence for neutrinos masses in this parameter range inside the PNS core. 

On the other hand, if one takes into account the KamLAND results 16 , which 
definitively demonstrated that i) a large mixing angle (LMA) solution of the solar v 
problem is favored: sin 2 26 <~ 0.8, ii) for a mass-squared difference: Am 2 ^ 5.5 x 10~ 5 
cV 2 (we use the approximate value Am 2 ~ 10~ 4 eV 2 for the estimates below), one 
can see that resonant conversions with Am 2 = m 2 2 — m 2 x <~ 10~ 4 eV 2 would 
take place in supernova regions where the density is about p <~ 10 — 30 g cm~ 3 , 
which corresponds to the outermost layers of the exploding star. Although a large 
number of v species can in effect participate of the transitions there, i.e, the neutrino 
luminosity can be still a very large quantity, these regions are of no interest for 
the gravitational-wave emission from neutrino oscillations since the overall energy 
density at that distance from the star center is very small. This does not mean 
that no GWs are emitted from transitions there, it is to mean that their strain is 
very small so as to be detectable. Observations of v oscillations 13 in this parameter 
range would provide useful information regarding the SNP, the hierarchy of neutrino 
masses and the mixing |C/ e 3| 2 . Note in passing that oscillations in this range would 
imply a mass for the v 2 species m„ 2 <~ 10~ 2 cV ( 10 g cm--' )> 111 tnc case wnen 
m„ 2 > m Vl . This to„ 2 is compatible with current limits from WMAP 20 . 

Finally, let us consider oscillations in the parameter range estimated from CMBR 
by WMAP observations. In this case, resonant v transitions would take place in 
regions where the density is as high as p <~ 10 7 ~ 8 g cm ~ 3 , that is, at the supernova 
mantle or PNS upper layers. At these densities the oscillation length can be still 
A osc ^1 — 5 km, and thus the conversions can be considered as adiabatic. Thus the 
resonance condition can be satisfied for Am 2 <~ 1 cV 2 as long as sin 2 26 m < 10~ 3 . 

Hence, by taking with this constraint, and recalling that i) at least 6 ^-species 
can participate in the flavor transitions, ii) most ^ s are emitted having k v parallel to 
B, which implies a further reduction factor of 2, and Hi) most vs are emitted having 
k v parallel to J implying also an additional reduction factor of 2, one can show from 
Eq.(3) that the fraction of v a species that can eventually exchange flavor during the 
first few milliseconds after core bounce turns out to be 3 P Va _ iVa (\x — x res \) < 0.1 % . 
Thus, the total energy involved in the oscillation process we are considering could 
be estimated as: £* ot ~ P„ a _> Va (\x - x, res \) x N v k v c (1 + 0.15), with N u the total 
number of neutrinos undergoing flavor conversions during the timescale AT osc . Note 
in passing that KamLAND D e experiment suggests 16 P Pe _>p M T <~ 40%, while LSND 
Pv^D e ~ 0.26%. This result from LSND has not been so far ruled out by any 
terrestrial experiment, and there is a large expectation that it could be verified by 
MiniBoone at Fermilab. 

In the case of active-to-sterile v conversions, it was proved in Refs. 3,4,1 (see 
details therein) that a substantial fraction: P Va _> Vs , of ^s may get converted to 
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sterile vs, and escape the core of the star, if the sterile v mass (m„ 4 ) is such that 

1 eV 2 < Am^„ 4 < 10 4 eV 2 . (6) 

For both classes of v conversions the number of ^s escaping and their angular 
distribution is sensitive to the instantaneous distribution of production sites. These 
inhomogeneities can give rise to quadrupole moments that generates GWs 3 ' 12 , and 
dipole moments that could drive the runaway pulsar kicks 5 ' 6 . Noting that at least 6 
^-species can participate in both types of oscillations and that the interaction with 
both the magnetic field and angular momentum of the PNS brings with an overall 
reduction factor of 4 in the oscillation probability, one can show that the fraction 
of vs that can actually undergo flavor transitions in the first few milliseconds is 3 

p^:(w-x \)<i%, (?) 

of the total us number: M x 1 ~ 10 57 ^, which corresponds to a total energy 
exchanged during the transition 

AE^ZZt. ~ 3 x 10 51 erg. (8) 

Equations (8,3) determine the total luminosity of the neutrinos participating in 
the resonant flavor conversions. 



4. GWs Energetics From u Luminosity and Detectability 

If v oscillations do take place in the SN core, then the most likely detectable GWs 
signal should be produced over the time interval for which the conditions for flavor 
conversions to occur are kept, i.e., ~ (10 _1 — 10) ms. This timescale implies GWs 
frequencies in the band: /gws ~ [10 - 0.1] kHz, centered at 1 kHz, because of the 
maximum v production around 1 — 3 ms after core bounce 3 . This frequency range 
includes the optimal ban width for detection by ground-based observatories. For a 
1 ms conversion time span the v luminosity reads 

3 x 10 51 erg 



AE V 
AT, 



osc 



1 x 10" 3 s 



= 3 x 10 



r,i 



erg 



(9) 



To estimate the GWs burst amplitude produced by the non-spherical outgoing 
front of oscillation-produced v s one can write (see details in Refs. 1 ' 3 ' 11 ' 12 ) 



2G 

hv = ~c T R ^ At X Lv X °^ ~ 



55kpc 



R 



AT 

IF 3 ! 



3 x 10 54 ^ 



\0.lJ 



where \A\ = 4 x 10 23 Hz 1/2 is the amplitude. b A GWs signal this strong will 
likely be detected by the first generation of GWs interferometers as LIGO, VIRGO, 



b One must realize at this point that the high value of the anisotropy parameter here used is 
consistently supported by the discussion by Mosquera Cuesta and Fiuza 1 regarding the neutrino 
coupling to rotation and magnetic field quoted above. 
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etc. Its imprint in the GWs waveform may resemble a spike of high amplitude and 
time width of <~ ms followed by a Christodoulou's memory 4 . From Eq.(10) the GWs 
luminosity turns out to be 

-i 2 / \ 2 



J GWs 



10 



,48-50 er S 



3 x 10 54 ^£ 



0.1 - 0.01 



(10) 



while the GWs energy radiated in the process yields ^GWs = ^GWs x ^^bsc ~ 

1Q 47-45 

erg . This is about 10 5 3 x L^ iS - the luminosity from v diffusion inside the 

PNS 4 . 



5. Conclusion 

One can see that if v flavor conversions indeed take place during SN core bounce 
inasmuch as they take place in our Sun and Earth 23 , then GWs should be released 
during the transition. The GWs signal from the process is expected to irradiate 
much more energy than current mechanisms figured out to drive the NS dynamics 
at birth do. A luminosity this large (Eq.(9)) would turn these bursts the strongest 
GWs signal to be detected from any SN that may come to occur on distances up 
to the VIRGO cluster, R ~ [10 - 20] Mpc. It is stressed that this signal will still be 
the stronger one from a given SN, even in the worst case in which the probability of 
v conversion is three orders of magnitude smaller than the estimated in the present 
paper. In proviso, we argue that a GWs signal that strong could have been detected 
during SN1987a from the Tarantula Nebula in the Milky Way's satellite galaxy 
Large Magellanic Cloud, despite of the low sensitivity of the detectors at the epoch. 
In such a case, the GWs burst must have been correlated in time with the earliest 
arriving neutrino burst constituted of some active species given off during the very 
early oscillation transient where some v e s went into v^s, v T s or ^ s s. Thenceforth, it 
could be of worth to reanalyze the data collected far from that event taking careful 
follow up of their arrival times, if appropriate timing was available at that moment. 
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